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The generalized Thomas-Kuhn sum rules are used to eliminate the explicit dependence on dipolar 
terms in the traditional sum-over-states (SOS) expression for the second hyperpolarizability to 
derive a new, yet equivalent, SOS expression. This new dipole-free expression may be better suited 
to study the second hyperpolarizability of non-dipolar systems such as quadrupolar, octupolar, 
and dodecapolar structures. The two expressions lead to the same fundamental limits of the off- 
resonance second hyperpolarizability; and when applied to a particle in a box and a clipped harmonic 
oscillator, have the same frequency-dependence. We propose that the new dipole-free equation, when 
used in conjunction with the standard SOS expression, can be used to develop a three-state model 
of the dispersion of the third-order susceptibility that can be applied to molecules in cases where 
normally many more states would have been required. Furthermore, a comparison between the two 
expressions can be used as a convergence test of molecular orbital calculations when applied to the 
second hyperpolarizability. 

PACS numbers: 42.65.An, 33.15.Kr, 11.55.Hx, 32.70.Cs 



I. INTRODUCTION 

The sum-over-states (SOS) expressions have been used for more than three decades in the study of nonlinear optical 
phenomena, and are perhaps the most universally used equations in molecular nonlinear optics. The sum-over-states 
expression is obtained from quantum perturbation theory and is usually expressed in terms of the matrix elements 
of the dipole operator, —exnm, and the zero-field energy eigenvalues, -En-fJ, @, Q 

The SOS expressions for the first and second hyperpolarizability derived by Orr and Ward using the method of 
averages are often used because they explicitly eliminate the unphysical secular terms that are present in other 
derivations. [ij These secular-free expressions contain summations over all excited states. 

Finite-state approximations are used to apply the theory to experimental results. Oudar and Chemla studied the 
first hyperpolarizability of nitroanilines by considering only two states, the ground and the dominant excited state. Q 
Although the general validity of this "two-level" model has been questioned, especially in its use for extrapolating 
measurement results to zero frequency, the approximation is still widely used in experimental studies of the nonlinear 
properties of organic molecules. 

Several approaches have been used to develop approximate expressions for the second-hyperpolarizability in 
the off-resonance regime. 0, 0, While such approximations are helpful, they systematically ignore some of the 
contributions to the SOS expression. As our goal is to derive a general expression that is equivalent to the traditional 
SOS one, we choose not to make any assumptions a priori about what type of contributions dominate the response. 
Furthermore, including all the possible contribution is necessary to properly describe the on-resonance behavior, even 
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when only few states contribute to the response. [8[ 

In 2005, Kuzyk used the generahzed Thomas-Kuhn sum rules to relate the matrix elements and energies involved 
in the general Orr and Ward SOS expression for the first hyperpolarizability, and introduced a new and compact 
SOS expression that does not depend explicitly on dipolar terms. ^] Since the Thomas-Kuhn sum rules are a direct 
and exact consequence of the Schrodinger equation when the Hamiltonian can be expressed as H = /2m + y(r), it 
follows that the new SOS expression is as general as the original, converges to the same results, and by virtue of its 
compactness may be more appropriate for the analysis of certain nonlinear optical properties. [lo| Indeed, Champagne 
and Kirtman used a comparison between the dipole-free and standard SOS expressions to study the convergence of 
molecular-orbital calculations. pD| In this work, we use the same principle to derive a compact and general dipole-free 
expression for the second hyperpolarizability. 



II. THEORY 



While our method can be applied to non-diagonal components of the second hyperpolarizability, for simplicity we 
will focus on the diagonal component. The SOS expression for the diagonal term of the second hyperpolarizability 7 
as derived by Orr and Ward in 1971 is given by:[2| 



' _ _ 00 ' 

where e is the magnitude of the electron charge, Xnm the n, m matrix element of the position operator and huji 
{i = 1, 2, 3) are the frequencies of the photons with lo^ =^^1+^2+^3. The bar operator is defined as: 
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The dispersion of 7 is given by /^^"^^(cji, 0^2, W3) and I?^^(cji, 1.^2, (j^s) which are defined as follows: 
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where spontaneous decay is introduced by defining complex energies: 

hfln = Eno - iVn, (5) 

where Eno is the energy different between the n*'' excited state and the ground state, and ^ is the inverse radiative 
lifetime of the n*'* state. 



A. Dipole-free expression for the second hyperpolarizability 



To obtain a dipole-free expression for the second hyperpolarizability we begin by separating explicitly dipolar terms 
from dipole-free terms in the first term of Eq. [1] 
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The second term in Eq. [l]is already dipole-free. 



(6) 



It should be noted that for non-dipolar systems (such as octupolar chromophores) , with Axmo = 0, only the last 
term in Eq. [S] contributes to the second hyperpolarizability. The generalized Thomas-Kuhn sum rules can be used to 
obtain a relationship between the explicitly dipolar terms in terms of only non-dipolar terms: Q 

I i2a sr^' {Enk + Eno) , , 

Ffcol AXkO = - V XokXknXnO- (7) 

■~: ^kO 

n^k 

We stress that the only assumption made in the derivation of Eq. [7] is that the sum rules hold, which is the case 
when the unperturbed Hamiltonian describing the system is conservative. 

Substituting Eq. [7] into Eqs. [6] and [1] yields the dipole-free expression for the second hyperpolarizability: 
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(8) 



So, Equation [8] is as general as the traditional sum-over-states expression. [1| 



III. APPLICATIONS 



It is useful to compare the convergence between the dipole-free expression for the second hyperpolarizability 
(Eq. [8|) with the traditional Orr and Ward SOS expression (Eq. [ij for various systems. In this section we will 
compare these expressions as a function of wavelength for two model systems. Mathematically, both expressions are 
equivalent, as long as all excited states of the system are included in the sum, so this exercise will determine how 
many states are required for convergence. Since in practice, the sum-over-states expressions must be truncated, it is 
critical to understand the effect of discarding terms on the nonlinear susceptibility. We also apply this new expression 
to calculate the fundamental hmits of 7, and show that the results agree with those obtained using the standard SOS 
expression. 
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A. Three-level model dipole-free expression: calculation of the fundamental limit in the off-resonance 

regime 

We begin by first calculating the fundamental limit of 7 starting from the dipole-free expression. The analogous 
calculation has already been performed using the traditional Orr and Ward SOS expression, so we can check 
whether or not the two results are the same. A different set of results would suggest that the method used in 
calculating the fundamental limits does not hold. 

According to the three-level ansatz,[l^[l^, when near the fundamental limit, only three-levels contribute to the 
nonlinear response, Eq. [8] becomes: 

lxxxxi-^a'-,^l,'^2,^3) f -Dill (2-^20 - -Eio)^ , ('^^20 ~ Eif)) \ ^ ^ 
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Off-resonance, the dispersion terms (Eqs. [3]and|4]) simplify to: 



= 4 \ jr-T^ > ' (10) 
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and 



^™^/ = 2 Tj^ri^ + i;^^ > (11) 

and the relationships between the first transition dipole moments can be evaluated from the Thomas-Kuhn sum 
rules: [1, [H, [H [M iH [11 [11 

ixo2p<f^[i4r^p-ixoin (12) 
ki2p< ,., ""^v j kr^p+KiH, (13) 



with 
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where iV is the number of electrons in the system. 



(14) 



Introducing the dimensionless quantities: 



E=^, (15) 
Koi I 

the off-resonance diagonal component of the second hyperpolarizability can be written as: 

Ixxxx = (enr (^^y ^^F,{E,X), (17) 

where Fj{E,X) is defined by: 

F^{E, X) = -5{E - lf{E +1){E'^ +E + 1)X'^ - 2{E^ - l)E^X^ - {E^ + E + 3)E^ + 4. (18) 
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The second hyperpolarizability scales as 7 oc iV^, the square of the number of delocalized electrons, and as 7 oc A^^^., 
the fifth power of the wavelength of maximum absorption. While in a three-level model, the expression for the first 
hyperpolarizability, which is analogous to Equation[T8l explicitly separates into a product of a function of the transition 
dipole moment xqi and excited state energies (i.e. Fp{E,X) — G{X)f{E))J\BJ\^ this is not possible for the case of 
the second hyperpolarizability. To optimize the second hyperpolarizability the function F-y{E, X) has to be optimized 
as function of the two parameters E and X. The behavior of the function F^{E, X) (as given by Eq. [T8|) as a function 
of the parameters X and E is shown in Fig. [TJ The function is maximized when both E (i.e. the second excited 
state energy level is far away from the first excited energy level) and X (i.e. the oscillator strength is concentrated 
in the second transition dipole moment, |a;o2|'^)- When the fimction is optimized we obtain the quantum limit: 

7;— = 4(e.r(^)'^. (19) 
This result agrees with the quantum limit obtained from the traditional Orr and Ward SOS expression, [l^ 




FIG. 1: F^{E,X) (as given by Eq. llSp spanning the full allowable range of the parameters X and E (defined in Eqs. 1151 and 
I16|) . The function is optimized when simultaneously _B — > and X 0. As X 1 and i5 ^ or _B ^ 1 the function becomes 
negative, corresponding to a negative value of the second hyperpolarizability. 



Thus, we can conclude that when only three levels contribute to the nonlinear response, the dipole-free and the 
traditional SOS expressions for the second hyperpolarizability - when simplified using the sum rules - become the 
same, leading to the same quantum limits. We should point out that the quantum limits are obtained by assuming 
that the response is dominated by the contributions of three overlapping states, an ansatz that has been extensively 
verified numerically using Monte Carlo methods [l^l as well as potential energy optimization. f2l| There are no 
assumptions about the symmetry properties of the states. However, this does not imply that symmetry plays no role 
in he optimization of the second hyperpolarizability. The symmetry properties of the system will determine whether 
or not the optimal distribution of excited energies and transition dipole moments can be achieved. Mathematically, 
symmetries will impose further constraints on the parameters, which will make 7 smaller. We note that the quantum 
limit is negative for the centrosymmetric system, and is one-quarter of the positive limit that is obtained for an 
asymmetric molecule. 



B. The particle in a box 



In this section, we test the convergence of the expressions in the case of two exactly solvable quantum mechanical 
systems: the "particle in a box" and the "clipped harmonic oscillator" . For simplicity, we will first perform our 
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calculations in the ofF-resonance regime. 

The unperturbed states that we will use for our calculation of the second hyperpolarizability are the solutions of 
the one-dimensional time-independent Schrodinger equation: 

The potential that characterizes the particle in a box is zero inside the box of length L and infinite, otherwise. The 
solutions are given by: 

,pb _ . ( + l)7rx 



with fc = 0, 1, 2, 3, . . . The corresponding energies are: 

where m is the mass of the particle (in this case the electron mass). These solutions are substituted into the 
expressions for the diagonal component of the second hyperpolarizability (Eqs. [T]and[5]) to study the convergence of 
both series as a function of number of excited levels included in the sum. 

In the off-resonance regime, Eq. [T] becomes: 
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and Eq. [5] is given by: 
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(24) 



The numerical evaluation of the diagonal component of the second hyperpolarizability as given by Eq. [23] or Eq. [24] 
respectively, is performed by dividing every contribution in the sum by the quantum limit (Eq. I19p . For the particle 
in a box, a general term in the sum can be rewritten as: 

Tn,T, T T n 3^ ^ ^ ^ ^ 

4 ■^m-^lm-^mn-^nO _ ^ off-max ap; Ulm ymnHna (0K\ 

IP IP 4 Ixxxx ' rPB tPB J^PB ' 

where we have used defined the following dimensionless functions: 

9^nn = / sm{{rn+l)y) ■ysm[{n + l)y)dy, (26) 
= n{n + 2). (27) 

The convergence of the two series in the ofF-resonance regime is shown in Fig. [2] When including the contribution 
of the first 10 states, the relative difference between traditional SOS expression for the second hyperpolarizability (Eq. 
[1]) and the dipole-free expression (Eq. [8|) is of the order of 5 x 10^''. With 50 states, the two expressions converge to 
the same value of: 

ott-PB 

J^^^ = -0.08936. (28) 
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Interestingly, the average value (also shown in Fig. [2]) is more accurate when as few as three states are included in 
the sum. This is because, for the particular case of the particle in a box, if not enough states are included in the 
sums, the traditional SOS expression tends to underestimate while the dipole-free expressions tends to overestimate 
7. In this case, with few excited states, using the average value yields a more reliable estimation of the second 
hyperpolarizab ility . This same result was found for the first hyperpolarizability^] and for studies used in modeling 
real molecules. [ill 
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FIG. 2: Convergence of the normalized SOS expressions as a function of the number of states in the sum for the particle in a 
box. Both the traditional SOS expression (Eq. [ij and the dipole-free expression (Eq. [8} converge quickly to the exact value 
given by Eq. [28l For completeness, we also plot the normalized average value of the two expressions, which converges more 
rapidly as a function of the number of states than the traditional and dipole-free expressions alone. The inset shows a mgnified 
view of the region indicated by the dash boxes. 



C. The clipped harmonic oscillator 



Another exactly solvable system is the clipped harmonic oscillator, whose potential energy function is given by: [2 

if a; < 0, 



V{x) 



00 

ifx>0 



where uj has dimensions of frequency. 
Introducing the dimcnsionless variable: 



the solutions are expressed as: 



^cho 



if ^ < 0, 
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where fffc(^) is the fc* order Hermite Polynomial, and the energies are given by: 



(29) 
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with fc = 0, 1,2,3, 
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For the clipped harmonic oscillator, a general term in the sum can be rewritten as: 

ottCHO nCHO „CHO nCHO nCHO 

^4 XQlXlmXmnXnQ _ Ixxxx , 9^21+1)9 (21+1) (2ni+l)9 (2ni+l){2n+l)9 (2n+l)l 

EioEjnoEno l-m-n 
where we have defined the following dimensionless function: 

/•OC 

g^^o = (2""in!)-i/2(2™-i,„!)-i/2 / H,,{x)xH^{x)dx. (34) 



The convergence of the two series in the off-resonance regime is shown in Fig. [31 In this case, it takes more terms 
to reach convergence than for the particle in a box. For 50 states, the relative difference between the traditional SOS 
expression (Eq. [1]) and the dipole-free expression (Eq. [8]) is of the order of 2 x 10~^. The traditional SOS expression 
converges faster to the final value: 

otf-CHO 

'^'^^'^^ = -0.00842. (35) 
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FIG. 3: Convergence of the normalized SOS expressions as a function of the number of states in the sum for the clipped 
harmonic oscillator. Both the traditional SOS expression (Eq. [TJ and the dipole-free expression (Eq. [SJ converge slowly to the 
exact value given by Eq. |35] For completeness, we also plot the normalized average value of the two expressions. The inset 
shows a magnified view of the region indicated by the dash boxes. 



D. Dispersion studies using 6 excited levels 



In this section, we investigate the convergence of the dispersion of the two forms of the second hyperpolariz- 
ability. In particular, we will treat two separate cases: two photon absorption (TPA) - which is related to the 
imaginary part of 7 and used in a broad range of applications such as photodynamic cancer therapies and 3D 
photolithography, jlj, [2^, [2^ and the optical Kerr effect (OKE) - which is related to the re al p art of 7 and 
widely used in characterizing materials with potential applications in all-optical switching. [13, [H, [2^ We begin by 
evaluating Eqs. [3] and [Hand including only the first 6 excited states. 

To get the typical qualitative behavior of real large chromophores, we choose Eiq = leV. All the linewidths are 
given by r„ = O.leV which is also typical for organic chromophores. We normalize the second hyperpolarizability by 
dividing by the off-resonance limit for £^10 — leV. 

Fig. [31 shows the dispersion predicted by the traditional sum-over-states expression and the dipole-free expression 
for the imaginary part of the second hyperpolarizability as s function of the incident photon energy. Note that 
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FIG. 4: The normalized imaginary part of the two photon absorption second hyperpolarizability as a function of the incident 
photon energy for a 7-level model for the particle in a box using the traditional sum-over-states and the dipole-free expressions. 
The inset shows a magnified view of the region indicated by the dashed box. Eio = leV. 



Im['y] is related to the two photon absorption cross-section. The agreement between the two expressions is ex- 
cellent everywhere with the exception of the third resonance (see inset in Fig. where the two differ by less than 20%. 

Next we compare the two expressions for the real part of the second hyperpolarizability - which is related to the 
optical Kerr effect - as a function of the fundamental photon energy for the particle in a box . The results are plotted 
in Fig. O In this case, for 6 excited states, the two expressions differ by as much as a factor of 2 in some regions 
(see insets). Although the qualitative behavior is the same, this type of discrepancy should be considered when 
experimental data is analyzed using a limited number of terms in the sum-over-states expressions. 
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FIG. 5: The normalized real part of the Kerr effect second hyperpolarizability as a function of the incident photon energy for 
a 7-level model for the particle in a box using the traditional sum over states and the dipole-free expressions. The insets show 
a magnified view of the region indicated by the dashed boxes. Eio — leV. 



Next we consider the dispersion predicted by the traditional sum-over-states expression and the dipole-free expres- 
sion for the imaginary part of the two photon absorption second hyperpolarizability as a function of the fundamental 
photon energy for the clipped harmonic oscillator. The results are plotted in Fig. [S] The agreement between the 
two expressions is good, differing only near the resonances. Finally, Fig. [7] shows the dispersion predicted by the 
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FIG. 6: The normalized imaginary part of the two photon absorption second hyperpolarizability as a function of the incident 
photon energy for a 7-level model for the clipped harmonic oscillator using the traditional sum-over-states and the dipole-free 
expressions. The inset shows a magnified view of the region indicated by the dashed boxes. Eio = leV. 

two expressions for the Kerr effect second hyperpolarizability as a function of the fundamental photon energy for the 
clipped harmonic oscillator. For 6 excited states, the two expressions differ substantially only in the vicinity of a 
resonance (see insets). 
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FIG. 7: The normalized real part of the Kerr effect second hyperpolarizability as a function of the incident photon energy for 
a 7-level model for the clipped harmonic oscillator using the traditional sum-over-states and the dipole-free expressions. The 
insets show a magnified view of the region indicated by the dashed boxes. Eio — leV . 



IV. CONVERGENCE OF THE DIPOLE-FREE EXPRESSION 

In this section we will study the convergence of the dipole-free series expression (Eq. [8|) as a function of the number 
of excited states included in the sum. The convergence is studied as a function of photon energy for two photon 
absorption and the Kerr effect, using the particle in a box and the clipped harmonic oscillator as model quantum 
systems. Again, in order to get the typical qualitative behavior of large real chromophores, we choose Eiq = leV, 
and all the linewidths are given by r„ = O.leV which is also typical for organic chromophores. Also, we normalize 
the second hyperpolarizability by dividing by the off-resonance limit for Eio = leV. 
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In all the cases it is found that after including about 15 excited states in the sum the expression converges even in 
the vicinity of a resonance. In order to get a better understanding of the convergence behavior of the expressions we 
will use the particle in a box model as a test for convergence, once we take into account symmetry considerations. 

For clarity, we will look again to the expression for the second hyperpolarizability that separates explicitly dipolar 
terms from dipole-free terms: 

oo'/oo'/oo' _ _ \\ 
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For a system whose symmetry demands that all dipole moments vanish, i.e. Axmo — for m — 1,2,3, ■• ■, such 
as the the centrosymmetric system of a particle in a box or molecules with octupolar symmetry (i.e. no dipole 
moment but non-centrosymmetric) , the first three terms must each vanish. Numerically, when we use the dipole-free 
expression to model such systems that are centrosymmetric, which demands that all dipole moments vanish, the 
contributions from the first three terms are precisely zero only when an infinite number of terms are included in the 
sum. 



It is useful to consider how many states are needed in order for each of these dipolar terms, written in our new 
non-dipolar form, to vanish. We define the partial sums 5*1, ^2 and S3 and ^4 as: 
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Clearly, for systems where there is no change in the dipole moment between the ground and excited 
states, the dipole-free expression will converge when simultaneously: 5*1 — > 0, S'2 — > and S3 and 
'S'4 — > jxxxx{—^a', i-^It '-^2i ^3) ■ We will test these conditions numerically using the particle in a box as a model of a 
centrosymmetric potential with no dipole moment. 

First we consider two photon absorption and study the convergence of Si (Eq. [27|), S2 (Eq. [35]), ^3 (Eq. [33]) and S'4 
(Eq. |40)) as a function of the number of excited states for a range of different photon energies. The results are shown 
in Fig. [5] Clearly, for all cases, few excited states (from 2 to 5) are needed to get good convergence when away from 
resonance. While Si also converges in the resonant regime after 10 excited states are included more excited states 
(up to 20) are required for convergence of ^2 and 6*3 near resonance (see insets). Finally we look at the convergence 
of S'4 (Eq. I40p . which as we have seen, converges to the exact value of the second hyperpolarizability for systems 
with no permanent dipole moment, such as the particle in a box. Surprisingly, this expression is shown to converge 
rapidly as a function of number of excited states included in the sum, even on resonance. In fact, it is clear from the 
plot that only 2 excited levels might suffice to study the qualitative behavior of the second hyperpolarizability even 
close to the resonances. 



Similarly, we next consider the Kerr effect and study the convergence of Si (Eq. [37|) , S2 (Eq. 135)) and S3 (Eq. 135)) 
and S4 (Eq. 1^0)) as a function of the number of excited states for a range of different photon energies. The results are 
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FIG. 8: The normalized imaginary part of the partial sums Si (top- left), S2 (top-right), S3 (bottom- left) and S4 (bottom-right) 
for two photon absorption as a function of the incident photon energy. The number of states included in the sum is shown 
in the figure legends. The particle in a box is used as a model centrosymmetric system with vanishing dipole moments, and 
therefore as the sums converge we must have Si 0, S2 ^ 0, S3 ^ and Im{S4,) = Iin{-f]:^^^^^). The insets show a magnified 
view of the regions indicated by the dashed boxes. 



shown in Fig. [51 Again, few excited levels are needed for convergence away from resonance, although more excited 
states (up to 20) are required for convergence near resonance. As in the case of the two photon absorption second 
hyperpolarizability, ^4 (Eq. I40p . which is the only terms that contributes to the second hyperpolarizability for the 
centrosymmetric potential given by the particle in a box is shown to converge rapidly as a function of number of 
excited states included in the sum, even in the resonant regime. 
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FIG. 9: The normalized real part of the partial sums Si (top- left), S2 (top-right), S3 (bottom- left) and S4 (bottom-right) for 
the Kerr efi'ect second hyperpolarizability as a function of the incident photon energy. The number of states included in the 
sum is shown in the figure legends. The particle in a box is used as a model centrosymmetric system with vanishing dipole 
moments, and therefore as the sums converge we must have Si ^ 0, S2 ^ 0, S3 — > and Re{S4,) = Rei'y^^^^'^^). The insets 
show a magnified view of the region indicated by the dashed boxes. 

For completeness, we also study the convergence of the expressions for the two photon absorption process and 
the Kerr effect using the clipped harmonic oscillator as a quantum model. As shown by the results in Fig. [TO] the 
expressions do not converge as rapidly as did the particle in a box; but when 15 excited states are included in the 
sum, the spectral features do not change quantitatively. 

All of our results show that the dipole-free expression converges even in the resonant regime when enough excited 
states (up to 20) are included in the expression. We have also shown that for systems with no dipole moments (such 
as octupolar molecules or the particle in a box) the dipole-free expression for the second hyperpolarizability collapses 
to Eq. [301 (which we will call the reduced dipole-free expression), since all the other terms vanish. In this case, the 
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FIG. 10: The normalized imaginary part of the two photon absorption second hyperpolarizability (top) and the normalized real 
part of the Kerr effect second hyperpolarizability(bottom) as a function of the incident photon energy. The number of states 
included in the sum is shown in the figure captions. The clipped harmonic oscillator is used as a model, with Eio = leV. The 
insets show a magnified view of the regions indicated by the dashed boxes. 

reduced dipole-free expression converges when only a few excited states are included in the sum. 

V. CONCLUSIONS 

We have developed an expression that eliminates the explicit dependence on dipolar terms but is physically 
equivalent to the traditional SOS expression for the second hyperpolarizability. The equivalence between the 
dipole-free and the traditional SOS expressions is demonstrated by calculating the quantum limits and studying the 
convergence of the series with the exact wavefunctions of two quantum systems: the particle in a box and the clipped 
harmonic oscillator. In both cases, when a large number of states is included, the two expressions are identical. 
However, the average of the two expressions converges faster than the individual expressions. 

Since the average between the two expressions appears to be a better approximation to molecular dispersion, the 
average may make it possible to use limited-state models when interpreting experimental dispersion data. Since 
accurate measurements of transition moments between excited states are difficult and tedious, the averaged second 
hyperpolarizability can be a useful tool for modeling the second hyperpolarizability when only limited information is 
available about the excited states of a particular system. 

To test the convergence between the two expressions, we have evaluated them in the resonant regime in two model 
systems: the particle in a box - which is a symmetric potential with no change in dipole moment; and the clipped 
harmonic oscillator - an asymmetric potential. This allows us to determine the role of symmetry. In both cases, 
we study the dispersion of the second hyperpolarizability near resonance, where - based on the different energy 
denominators - one would expect the differences between the two expressions to be the least consistent. The reduced 
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dipole-free expression has been introduced for systems with no dipole moment. Such an expression might be most 
appropriate when experimental results are interpreted since it requires the minimum number of molecular parameters. 

In conclusion, the dipole-free expression is an alternative to the traditional SOS expression that increases the 
theoretical pallet available to quantum chemists. It is more direct in certain theoretical problems such as its 
application to the derivation of a more-rigorous calculation of the fundamental limits of the third-order susceptibility. 
It provides a tool to assess the convergence of truncated SOS calculations, can be used to determine the accuracy of 
molecular-orbital calculations of nonlinear susceptibilities, and can be used to refine limited-state models to interpret 
experimental results. And, it may be more naturally applicable to the analysis of specific systems such as octupolar 
structures. [M ill [11 il [M, [11, [13, il 
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